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Abstract. In this note we announce that under general hypotheses, wavelet- 
type expansions (of functions in L p , 1 < p < oo, in one or more dimensions) 
C*| ■ converge pointwise almost everywhere, and identify the Lebesgue set of a func- 

tion as a set of full measure on which they converge. It is shown that unlike 
the Fourier summation kernel, wavelet summation kernels Pj are bounded by 
radial decreasing L 1 convolution kernels. As a corollary it follows that best 
L 2 spline approximations on uniform meshes converge pointwise almost every- 
where. Moreover, summation of wavelet expansions is partially insensitive to 
order of summation. 

We also give necessary and sufficient conditions for given rates of conver- 
gence of wavelet expansions in the sup norm. Such expansions have order of 
convergence s if and only if the basic wavelet ip is in the homogeneous Sobolcv 
space Hu S ■ We also present equivalent necessary and sufficient condi- 



tions on the scaling function. The above results hold in one and in multiple 
dimensions. 



1. Introduction 



In this note we present several convergence results for wavelet and multiresol- 
ution-type expansions. It is natural to ask where such expansions converge (and 
whether they converge almost everywhere) and what are the rates of convergence. 
The answer is that under rather weak hypotheses, one- and multidimensional 
wavelet expansions converge pointwise almost everywhere and, more specifically, 
on the Lebesgue set of a function being expanded. We will also give exact rates 
of convergence in the supremum norm, in terms of Sobolev properties of the basic 
wavelet or of the scaling function. 
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Wavelets with local support in the time and frequency domains were defined by 
A. Grossman and J. Morlet [GM] in 1984 in order to analyze seismic data. The 
prototypes of wavelets, however, can be found in the work of A. Haar [Ha] and the 
modified Franklin systems of J.-O. Stromberg |StJ. 

To identify the underlying structure and to generate interesting examples of 
orthonormal bases for i 2 (R), S. Mallat |Ma| and Y. Meyer [Me| developed the 
framework of multiresolution analysis. P. G. Lemarie and Y. Meyer |LM| con- 
structed wavelets in <S(R n ), the space of rapidly decreasing smooth functions. J.- 
O. Stromberg ptj developed spline wavelets while looking for unconditional bases 
for Hardy spaces. G. Battle [Ba] and P. G. Lemarie [Le| developed these bases 



in the context of wavelet theory. Spline wavelets have exponential decay but only 
C N smoothness (for a finite N depending on the order of the associated splines). 
I. Daubechies [Dal] constructed compactly supported wavelets with C N smooth- 
ness. The support of these wavelets increased with the smoothness; in general, to 
have C°° sm ooth ness, wavelets must have infinite support. 

Y. Meyer [Me] was among the first t o study con vergence results for wavelet ex- 
pansions; he was followed by G. Walter [Wal, Wa2 |. Meyer proved that under some 
regularity assumptions on the wavelets, wavelet expansions of continuous functions 
converge everywhere. In contrast to these results, the pointwise convergence re- 
sults presented here give almost everywhere convergence (and convergence on the 
Lebesgue set) for expansions of general LP (1 < p < oo) functions. We assume 
rather mild bounds and no differentiability for the wavelet or the scaling function; 
our conditions allow inclusion of the families of so-called r-regular wavelets [Mc|, 
as well as some other wavelets. 

These results parallel L. Carleson's 



Ca| and R. A. Hunt's [Hu] theorems for 



Fourier series. One difference and slight advantage of wavelet expansions comes 
from the fact that almost everywhere convergence occurs on a simple set of full 
measure, namely the Lebesgue set, while almost everywhere convergence for Fourier 
series is established on a much more elaborate set of full measure. 

We also give necessary and sufficient conditions for given pointwise (sup-norm) 
rates of convergence of wavelet or multiresolution expansions, in terms of Sobolev 
conditions on the basic wavelet or the scaling function. It has been shown previously 
by Mallat |Ma[ ] and Meyer [Me that the Sobolev class of a function is determined 
by the L 2 rates of convergence of its wavelet expansion. Necessary and sufficient 
conditions for L? rates of convergence which are analogous to our sup-norm con- 



ditions have been obtained by de Boor, DeVore, and RonjBDR], who have also 
studied sup-norm convergence in more general situations [BR . 

Our results on convergence rates ca n be viewed as a sharpening in the context of 
wavelets of the well-known Strang-Fix [pF| conditions for convergence of multiscale 
expansions. 

The results given here hold for multiresolution expansions in multiple dimension. 
The proofs, which will appear elsewhere [ KKR |, involve the kernels of the partial 
sums of such expansions and the above-mentioned result that such kernels are 
bounded by rescalings of L 1 radial functions. We should add that such bounds for 
wavelet expansions are nontrivial and arise from cancellations which occur in the 
sum representations of the partial sum kernels. Naive bounding of the summation 
kernels by using absolute values in their sum representations fails to yield the 
needed radial bounds for any class of wavelets. We remark that such bounds on the 
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summation kernel can be obtained more easily by writing it using the orthonormal 
translates <p(x — k) of the scaling function, instead of the wavelets ipjk- However, in 
proving results for convergence of wavelet expansions (Theorem 2.1 (hi)), we wish 
to avoid making any assumptions about radial bounds for the scaling function. 

The pointwise and L p convergence results contained here were obtained inde- 
pendently by the first author and by the second two authors. Results on the Gibbs 
effect obtained by the first author will appear elsewhere. 



To start, we define a multiresolution analysis on L (R ) Ma Vic] 



Definition 1.1. A multiresolution analysis on L 2 (R d ), d > 1, is an increasing 
sequence {Vj}j^z, 

• • • C V-2 C V-i C V C Vi C V 2 C • • • 

of closed subspaces in L 2 (R d ) where 



and the spaces Vj satisfy the following additional properties: 
(i) For all / S L 2 (R d ), j € Z, and k e Z d , 

f{x) e Vj <=> f(2x) e V j+1 

and 

f{x) £7»^ fix -k)eV . 

(ii) There exists a scaling function <p s Vq such that {4>jk}k£Z d is an orthonormal 
basis of Vj , where 

(1) ^ k {x) = V d ' 2 4>(^x-k), 

for x e R d , j G Z, and k E Z d . 

Associated with the Vj spaces, we additionally define Wj to be the orthogonal 
complement of Vj in V J+1 , so that V j+1 = Vj © Wj. Thus, L 2 (R d ) = J2® w j- We 
define Pj and Qj = Pj+i — Pj, respectively, to be the orthogonal projections onto 
the spaces Vj and Wj, with kernels Pjix,y) and Qjix,y). 

Under the assumptions in the above definition and with some additional regu- 
larity, it can be proved [Mc, Da2| that there then exists a set {ip x } & Wo, where 



A belongs to an index set A of cardinality 2 d — 1, such that {4 ! fk}kez d .x is an 
orthonormal basis of Wj, and thus {^jk}jez,kez d ,\ is a wavelet basis of L 2 (R, d ), 
where 

(2) Vi(z) = 2 J ' d/ V(2 J z-*0, 

for x e R d , j E Z, and k e Z d . 

Definition 1.2. For / G L p (R, d ) (1 < p < co) we define the following related 
expansions: 

(a) The sequence of projections {Pjfix)}j will be called the multiresolution 
expansion of /. 
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(b) The scaling expansion of / is defined to be 

(3) / ~ X] bk ^ x ~ fc ) + X! a J^jk( x ) ! 

fe i>0;fc:A 

where the coefficients Ojfc and && are the L 2 expansion coefficients of /. 

(c) The wavelet expansion of / is 

(4) /~E s 4w. 

i;fc;A 

where the coefficients aj/, are the L 2 expansion coefficients of /. 

We remark for part (a) of the above definition that it can be shown that the 
projections Pj (defined by their integral kernels) extend to bounded operators on 
L p , 1 < p < oo. The L 2 expansion coefficients in (b) and (c) (defined by integration 
against /) arc defined and uniformly bounded for any / G L p , 1 < p < oo. 

Definition 1.3. The point x is a Lebesgue point of the measurable function /(x) 
on R d if / is integrable in some neighborhood of x, and 

where B e denotes the ball of radius e about the origin, and V denotes volume. 

Such points p are essentially characterized by the fact that the average values of 
f(x) around p converge to the values of / at these points, as averages are taken over 
smaller balls centered at x. Note that all continuity points are Lebesgue points, 
but the converse is not true. 

Definition 1.4. A function f(x) is in the class 1ZB if it is absolutely bounded by 
an L 1 radial decreasing function r)(x), i.e., r)(xi) — n(x2) whenever xi = |x2|, 
v( x i) < v( x 2) when |xi| > |x2|, and r)(x) G L 1 (R d ). 

A function / is partially continuous if there exists a set A of vectors a G R d with 
positive measure such that lim e ^o ip{x + ea) = ifi(x) for a G A. 

Definition 1.5. The homogeneous Sobolev space of order s is defined by 



(5) H> h = Ue L 2 (R d ): ||/|| M ee J J \f(On\ 2s dC < oo \. 

The ordinary Sobolev space H s is defined as H^ in (5), with replacement of \£\ 2s 
by (|^| 2s + 1). Under the Fourier transform, the space H^ is a dense subspace of the 
complete weighted L 2 space of all measurable /(£) with ||/||/i, s < oo. This dense 
subspace consists of those functions /(£) which are also in the regular unweighted 
L 2 space. 

Convergence rates of wavelet expansions are sensitive to both the smoothness 
of the wavelet and the smoothness of the function being expanded. For a wavelet 
ip of given smoothness, the sensitivity to the Sobolev space of the function being 
expanded disappears when the function's Sobolev parameter is sufficiently large. 



POINTWISE CONVERGENCE OF WAVELET EXPANSIONS 5 

Definition 1.6. A family -0 A of wavelets yields pointwise order of approximation 
(or pointwise order of convergence) r in the space H s if for any function / G H s , 
the jth order wavelet approximation Pjf satisfies 

(6) ||P j /-/|| 00 =0(2-n 

as j tends to infinity. More generally, the wavelets ip yield pointwise order of 
approximation (or convergence) r if for any function / which is sufficiently smooth 
(i.e., is in a Sobolev space of sufficiently large order s), (6) holds. 

We will give several necessary and sufficient conditions (in terms of their Fourier 
transforms and membership in homogeneous Sobolev spaces) on the basic wavelet 
%j) or the scaling function <j> for given orders of convergence. In practice, sufficient 
smoothness for a function / (in the sense of the above definition) will mean that / 
is in H s + d / 2 or a higher Sobolev space. 

2. Pointwise convergence results 

With the background given, our main results can now be presented. 

Theorem 2.1. (i) Assume only that the scaling function cj> of a given multiresolu- 
tion analysis is in 1ZB, i.e., that it is bounded by an L 1 radial decreasing function. 
Then for any f G L p (H d ) (1 < p < oo), its multiresolution expansion {Pjf} con- 
verges to f pointwise almost everywhere. 

(ii) If (j>,ip G HB for all A, then also both the scaling (3) (if 1 < p < oo) and 
wavelet (4) (if 1 < p < oo) expansions of any f G L p converge to f pointwise almost 
everywhere. If further ip x and <fi are (partially) continuous, then both expansions 
converge to f on its Lebesgue set. 

(iii) If we assume only -0 A (x)ln(2 + |x|) G TVS for all X, then for f G L p , its 
wavelet (for 1 < p < oo) and multiresolution (for 1 < p < oo) expansions converge 
to f pointwise almost everywhere. If further ip x is (partially) continuous for all A, 
then both the wavelet and multiresolution expansions converge to f on its Lebesgue 
set. 

(iv) The last two statements hold for any order of summation in which the range 
of the values of j for which the sum over k and X is partially complete always 
remains bounded. 

In statement (iv) above, the summation over k and A is partially complete for a 
fixed j if it contains some terms, but not all with the given value of j. By the range 
of values for which the sum is partially complete we mean the difference of the 
largest and smallest value of j for which the sum is partially complete. Statement 
(iv) requires that this range always be smaller than some constant M. 

The above result on convergence of multiresolution expansions applies to spline 
expansions as well. Given a uniform grid K in R, one might ask whether given a 
function / G L 2 (R d ), the best L 2 approximations Pjf of / (by splines of a fixed 
polynomial order k) converge to / pointwise as the grid size goes to 0. The answer 
to this is affirmative. 

Corollary 2.2. For f G L P (R) (1 < p < oo), the order k best I? spline approxi- 
mations Pj f of f converge to f pointwise almost everywhere, and more specifically 
on the Lebesgue set of f, as the uniform mesh size goes to 0. 
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The proof follows from the fact that best spline approximations are partial sums 
of multiresolution expansions, with some radially bounded (1ZB) scaling function 
(f>. This result also extends to multidimensional splines. Technically, the best L 2 
approximation of / £ L p only makes sense when p = 2, but it can be defined for 
functions in L p by continuous extension of the projections Pj from L 2 to L p . 

It 



The following proposition is a consequence of the proof of Theorem 2.1 



has been proved before under somewhat stronger hypotheses, yielding stronger 



conclusions in [Me 



Proposition 2.3. Under the hypotheses of case (i), case (ii), or case (iii) of The- 
orem [2.l| ; L p convergence of the expansions also follows for 1 < p < oo. 

Thus for wavelet series and more generally for one- and multidimensional mul- 
tiresolution expansions, essentially all hoped for convergence properties hold, re- 
gardless of rates of convergence. 

The basis for Theorem 



2.1 



is the bound on the kernel of the projection Pj onto 
the scaling space Vj . It can be shown that under any of the hypotheses in Theorem 
2.1, the kernel P m (x,y) has the form 



P m {x,y) = Y^ ^jk&Wjkiy) = ^2 4>™k{x)<t>mk{y) 

j<m;k£Z d ;\ feGZ 

for x,y £ R d , with convergence of both sums on the right occurring pointwisc, 
uniformly on subsets a positive distance away from the diagonal D = {(x, y) : x — 
y}. The kernel converges to a delta distribution S(x — y) in the following sense: 

Theorem 2.4. Under the assumption that <fi £ 1ZB or that ip x (x) ln(2 + |x|) £ 1ZB 
for all X, the kernels Pj{x,y) of the projections onto Vj satisfy the convolution 
bound 

(7) \P 3 (x iy )\<Cy d H(y(\x~y\)), 

where H(\ ■ |) G 1ZB, i.e., H(\ ■ |) is a radial decreasing L 1 function. 

In one dimension, precise bounds can be obtained for kernels of specific wavelets. 
Two examples are illustrated in the result below. 

Theorem 2.5. InR 1 , let Pj(x,y) — X^cz 4>jk(x)4>jk(y) be the summation kernel, 
generated by the scaling function <fi £ L 2 (R). 

(a) If <j) has exponential decay, i.e., 4>{x) < Ce~ a ' x ' for some positive a, then 

\Pj(x,y)\ <CVe" a2l \ x - y \' 2 . 

(b) If (j) has algebraic decay, <fi(x) < ,. fj",^ for some N > 1, then 



\^y)\<c N {1 + 2 *_ y])N <c N 2> 



forN > 1. 
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As a corollary, if a scaling function (p has rapid decay (faster than any polyno- 
mial), then Pj(x,y) is bounded by a scaled convolution kernel (of the form on the 
right side of (7)) which has rapid decay. Spline wavelets (see |St], Ba, Lq]) satisfy 



the conditions in part (a), and wavelets constructed by P. G. Lemarie and Y. Meyer 
|LM| are an example of wavelets of rapid decay. 

3. Rates of convergence 

We will now give necessary and sufficient conditions on the basic wavelet ip x 
and on the scaling function </>, for given supremum-norm rates of convergence of 
wavelet expansions. The conditions on xjj (here given in terms of membership 
in homogeneous Sobolev spaces) can be translated into differentiability and then 
moment conditions on ip x (see the introduction). 

Theorem 3.1. Given a multiresolution analysis with either (i) a scaling function 
4> G 1ZB, (ii) basic wavelets satisfying if} X ln(2 + |ir|) 6 1ZB for each X, or (iii) a 
kernel for the basic projection P satisfying \P(x,y)\ < H(\x — y\) with H G 1ZB, 
the following conditions ((a) to (e')) are equivalent: 

(a) The multiresolution expansion (see Definition 1.2) yields pointwise order of 
approximation s. 

(a') The multiresolution expansion yields pointwise order of approximation s in 
every Sobolev space H r for r > s + d/2. 

(b) The projection I — Pj : H^ — * L°° is a bounded operator, where I is the 
identity and d denotes dimension. 

If there exists a family {i/j x }\ C 1ZB of basic wavelets corresponding to {Pj}, 
then: 

(c) For every family of basic wavelets {ip x }\ C 1ZB corresponding to {Pj}, and 
for each A, -0 A G H h s ' . 

(d) For every family of basic wavelets {ip*}\ C 7ZB corresponding to {Pj}, and 
for each X, 

(8) / i^(orir 2(s+d/2) ^<oo 

J\t\<e 

for some (or for all) e > 0. 

(d') For some family of basic wavelets {ip x }\ C 7ZB corresponding to {Pj}, 
equation (8) holds. 

If in addition there exists a scaling function (j) £ 1ZB corresponding to some 
family of basic wavelets as above, then : 

(e) For every such scaling function, 

(9) / ((2rr) d |^)| 2 -l)|er 2(s+d/2) ^<oo 

for some (or all) e > 0. 

(e') For some such scaling function <fi, equation (9) holds. 
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The above conditions on the scaling function can also be given in modified form in 
the case where the scaling function <f> has nonorthonormal integer translates. Neces- 
sary and sufficient conditions for convergence rates of spline and other nonorthonor- 
mal expansions can then be obtained directly from the same arguments as above. 
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